Testing General Relativity on Horizon Scales and the Primordial non-Gaussianity 
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The proper general relativistic description of the observed galaxy power spectrum is substantially different 
from the standard Newtonian description on large scales, providing a unique opportunity to test general relativity 
on horizon scales. Using the Einstein equations, the general relativistic effects can be classified as two new 
terms that represent the velocity and the gravitational potential, coupling to the time evolution of galaxy number 
density and Hubble parameter. Compared to the dominant density and velocity redshift-space distortion terms, 
the former scales d& T-L/k and correlates the real and imaginary parts of the Fourier modes, while the latter 
scales as {T-i/kY , where k is the comoving wave number and T-i is the conformal Hubble parameter. We use 
the recently developed methods to reduce the sampling variance and shot noise to show that in an all sky galaxy 
redshift survey at low redshift the velocity term can be measured at IO-ct confidence level, if one can utilize halos 
of mass M > lO^'^ Mq, while the gravitational potential term itself can only be marginally detected. We 
also demonstrate that the general relativistic effect is not degenerate with the primordial non-Gaussian signature 
in galaxy bias, and the ability to detect the primordial non-Gaussianity is little compromised. 
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The recent discovery of the cosmic acceleration has re- 
newed interest in modifications of gravity on cosmological 
scales, and tests of general relativity on horizon scales become 
ever more crucial in determining whether the cosmic acceler- 
ation is due to dark energy or the breakdown of general rela- 
tivity. However, there exists a fundamental limitation to any 
attempts for testing general relativity on large scales, and de- 
cisive conclusions remain elusive due to the cosmic variance 
limit set by the finite number of measurements available to us. 

In the past few decades galaxy clustering has been one 
of the indispensable tools in cosmology, covering a progres- 
sively larger fraction of the sky with increasing redshift depth. 
With the upcoming dark energy surveys this trend will con- 
tinue in the future. However, despite the advance in observa- 
tional frontiers, there remained a few unanswered questions 
in the theoretical description of galaxy clustering. On horizon 
scales, the standard Newtonian description naturally breaks 
down, and a choice of hypersurface of simultaneity becomes 
an inevitable issue, demanding a fully relativistic treatment of 
galaxy clustering beyond the current Newtonian description. 

In recent work 11101, it is shown that a proper relativistic 
description can be easily obtained by following the observa- 
tional procedure in constructing the galaxy fluctuation field 
and its statistics: We need to model observable quantities, 
rather than theoretically convenient but unobservable quan- 
tities, usually adopted in the standard method. While both the 
relativistic and the standard Newtonian descriptions are vir- 
tually indistinguishable in the Newtonian limit, they are sub- 
stantially different on horizon scales, rendering galaxy clus- 
tering measurements a potential probe of general relativity. 
The relativistic description of galaxy clustering includes two 
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new terms that scale as velocity and gravitational potential. 
Compared to the dominant density contribution, they are sup- 
pressed hy H/k and {T-L/kY and become important only on 
large scales, where the comoving wavevector amplitude k is 
comparable to the conformal Hubble parameter T-L ~ aH. 
Consequently, the identification of these terms just by looking 
at the galaxy power spectrum is hampered because of sam- 
pling variance, and the general relativistic effects unaccounted 
in the standard Newtonian description may result in system- 
atic errors less than l-a for most of the volume available at 
z < 3 in the standard power spectrum analysis 

A new multi-tracer method lU takes advantage of the fact 
that differently biased galaxies trace the same underlying mat- 
ter distribution, and it can be used to cancel the randomness of 
the matter distribution in a single realization of the Universe, 
eliminating the sampling variance limitation. This method has 
been used in ^ to investigate the velocity effects of 
noting that for any given Fourier mode the imaginary part of 
velocity couples to the real part of density and vice versa. 
Even with this novel technique, the expected detection level 
is low 0]. 

If sampling variance is eliminated, the dominant remain- 
ing source of error is shot noise, caused by the discrete nature 
of galaxies. Recently, a shot noise cancelling technique has 
been proposed IStl and investigated for detecting primordial 
non-Gaussianity |6!] in combination with the sampling vari- 
ance cancelling technique. The basis of the method is that 
by using halo mass dependent weights one can approximate 
a halo field as the dark matter field and reduce the stochas- 
ticity between them. While this works best when comparing 
halos to dark matter, some shot noise cancelling can also be 
achieved by comparing halos to each other This opens 

up a new opportunity to probe horizon scales and extract cos- 
mological information with higher signal-to-noise ratio. In 
this Letter we explore the possibility of using galaxy power 
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spectrum measurements combined with the muhi-tracer and 
shot noise cancelHng techniques to test general relativity on 
cosmological horizon scales. In addition, we comment on the 
impact of the general relativistic effects in detecting the pri- 
mordial non-Gaussian signature in galaxy bias. 

GR description. — A full general relativistic description of 
galaxy clustering is developed in |[ll|2] (see also iHIIId). Here 
we make one modification in the adopted linear bias ansatz. 
Previously, we have adopted the simplest linear bias ansatz, 
in which the galaxy number density is just a function of the 
matter density Ug = J^[pm], both at the same spacetime ifToll . 
However, this ansatz turns out to be rather restrictive, since 
the time evolution of the galaxy sample is entirely driven by 
the evolution of the matter density oc (1 + z)'^. Now we re- 
lax this assumption and provide more freedom, while keep- 
ing the locality: We allow the galaxy number density at the 
same matter density to differ depending on its local history, 
i.e., Ug = F[p„i ,tp] with tp being the proper time measured 
in the galaxy rest frame. Therefore, when expressed at the 
observed redshift z, the galaxy number density can be written 



ng{z) [l + b {Sm -3Sz) ~ ht Szy] 



ng{z) 



1 + 6 - e 5z^ 



(1) 



where the lapse 5z in the observed redshift z is defined as 
1 + z = (1 + 5z) I a and the expansion factor is a. The script v 
indicates quantities are evaluated in the dark matter comoving 
gauge (w = 0), and this gauge condition coincides with the 
synchronous gauge, as presureless dark matter can freely fall 
in the rest frame (see, e.g., lfl2ll ). 

The galaxy bias factor is 6 = d\Tingld\w.pra\t^ and the 
time evolution of the galaxy number density due to its local 
history is 6f = 91nng/91n(l + z)|p^. Therefore, the total 
time evolution of the mean gala xy number density is propor- 
tional to the evolution bias 11311 . e = o?lnng/(iln(l + z) = 



3 6 + ht. This biasing scheme in Eq. ([T]i is consistent with 
and our previous bias ansatz corresponds to 
ht = 0. Physically, the presence of long wavelength modes 
affects the local dynamics of galaxy formation by changing 
the local curvature and thus the expansion rate, and these are 
modulated by the Laplacian of the comoving curvature fl4ll . 

Therefore, with this more physically motivated bias ansatz, 
the general relativistic description of the observed galaxy fluc- 
tuation is flU 
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where the luminosity function slope of the source galaxy pop- 
ulation is p, the comoving line-of-sight distance is r, the fluc- 
tuation in the luminosity distance is SVl, the temporal and 
spatial metric perturbations are and Lp^, the line-of-sight 
velocity is V, and the gauge-invariant radial displacement and 
lensing convergence are STZ and /C. We have ignored the vec- 
tor and tensor contributions to (5gr ■ The subscript x indicates 
quantities are evaluated in the conformal Newtonian gauge 
(x = 0). We emphasize that compared to (J [2[| it is only 
the terms in the square bracket that are affected by the choice 
of linear bias ansatz. 

Noting that Eq. (|2|i can be most conveniently computed in 
the conformal Newtonian gauge, we first remove the lapse 
term Szu by using its gauge transformation property (5z„ = 
6z^ + Tivy^/k, but we keep the comoving gauge matter fluc- 
tuation to be consistent with the convention in the liter- 
ature. In Fourier space, a further simplification can be made 
by ignoring the projected quantities such as the gravitational 
lensing and the integrated Sachs-Wolfe contributions, which 
are important only for the pure transverse modes (/s" = 0) 
10, [i^ ■ Equation (|2]i in Fourier space is then 
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where / is the logarithmic growth rate and we defined the 
two dimensionless coefficients TZ and V in the last equality. 
We have used the Einstein equations in the second equality 
to express a^, ip^, and in terms of Sm^ . The standard 
Newtonian description of the observed galaxy fluctuation is 

S^cwt = b Sm^ - fil kv^/H = {b + fpl) Sm\ which can be 
contrasted with (5gr in Eq. (|2|. 



Apparent from their spatial dependence in Eq. (O, the coef- 
ficients TZ and V originate from the velocity and gravitational 
potential. While P is purely relativistic, some contributions 
to TZ may be considered non-relativistic, since they could be 
written down in Newtonian dynamics, simply as a coupling 
of velocity from the Doppler effect with the time evolution 
of the galaxy number density. However, correct estimates of 
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FIG. 1 ; Predicted measurement significance of general relativistic ef- 
fects in the galaxy power spectrum. All the halos above a given min- 
imum mass are utilized to take advantage of the multi-tracer method 
0,0]. Three curves represent different survey redshift ranges with 
corresponding volume V = 59, 410, 600 (/i^^Gpc)'^. Gray curves 
show the degradation of the measurement significance, when the evo- 
lution bias e and the luminosity function slope p are marginalized 
over. For the volume-limited sample with constant comoving num- 
ber density, (V, TV) = (1.7, 1.6), (2.1, 1.5), and (2.2, 1.5) at each 
redshift. 



the TZ value (and V) require fully relativistic treatments, as 
perturbations like velocity and potential are gauge dependent 
quantities. Furthermore, while Eq. Q can be derived with 
the minimal assumption that the spacetime is described by a 
perturbed FLRW metric and photons follow geodesies, the co- 
efficients TL and V in Eq. (O are obtained by utilizing the Ein- 
stein equations. Therefore, measuring TZ and V is equivalent 
to a direct measurement of the relativistic contributions, not 
through their influence on the matter density fluctuation, and 
we collectively refer to the coefficients TZ and V as the general 
relativistic effects in the galaxy power spectrum. 

Multi-tracer shot noise cancelling method.— We consider 
multiple galaxy samples with different bias factors for mea- 
suring the general relativistic effects in the galaxy power spec- 
trum. Using a vector notation, the observed galaxy fluctuation 
fields can be written as 



cobs _ 
"GR — 



bo + fi4i 



c-pV 

WW 



k/n 



(4) 



where bo, 7, e are the linear bias, the multi-dimensional iden- 
tity, and the residual-noise field vectors. Similarly, the coef- 
ficients TZ and 7-* become vectors with different values of e 
and p for each sample. We introduced two new parameters cji 
and c-p to generalize the measurement significance of the co- 
efficients TZ and 7^ to the case of multiple galaxy samples — 
they sue c-jz — c-p = 1 in general relativity, and measurements 
of these two parameters amount to the measurement signifi- 
cance of the two vectors TZ and 7-". By definition the noise- 
field is independent of the matter fluctuation {eSm'^ ) ~ 0, and 
the square bracket in Eq. (|4]i defines the effective bias vector b. 



In order to assess our ability to measure the general rela- 
tivistic effects in the galaxy power spectrum, we employ the 
Fisher information matrix: 



1 



-Re(7,y) 



Re - ap.^l 



(5) 



where a = b^£-^bP^, p, = b^£-^b,P„„ 7,, = 
blS^^bjPm, the shot noise matrix £ = (ee^), the matter 
power spectrum in the comoving gauge is Pm{k), the param- 
eters di = CTz,cp,e, p, and bj is the derivative of b with re- 
spect to the parameters 9i. This formula is a straightforward 
extension of the Fisher matrix derived in fdi], with the effec- 
tive bias vector b being a complex vector The imaginary part 
arises solely from the 7?.-term in Eq. (|4| and its derivative. 

Measurement significance. — For definiteness, we consider 
full sky surveys with three different redshift ranges and adopt 
a set of cosmological parameters consistent with the WMAP7 
results llTil . Given the survey volume V, the Fisher ma- 
trix in Eq. (|5]l is summed over the Fourier volume, where 
fcniin = 27r/V^i/3 and fc^ax = 0.03/iMpc"\ To model the 
Fisher matrix parameters a, 7^ , we adopt the halo model 
description in 10,01; It has been well tested against a suite of 
A^-body simulations with Gaussian and non-Gaussian initial 
conditions. We assume that the galaxy samples have a con- 
stant comoving number density (e = 37) and they are selected 
in a volume-limited survey to eliminate the magnification bias 
(or p = 0.47 in practice), but we marginalize over e and p in 
constraining c-n and c-p with priors de =0.1 and cTp = 0.05 
id. 

Figure [T| shows the predicted measurement significance of 
the general relativistic effects. In our most optimistic scenario, 
the velocity term c-ji can be measured at more than lO-cr con- 
fidence level at z < 1, while the gravitational potential cp is 
at 1-cr significance. At higher redshift, though the increase in 
the survey volume is partially cancelled by the lower abun- 
dance of halos at a fixed mass, a substantial improvement 
can be achieved by going beyond z — 1. However, as the 
gray curves show, the uncertainties in e and p need to be con- 
trolled beyond the current observational limit, before further 
improvement can be realized. 

Compared to the estimates S/N< 1 obtained by using the 
standard power spectrum analysis yj, the multi-tracer method 
dramatically enhances the measurement significance. Further- 
more, a method of measuring the imaginary part in the galaxy 
power spectrum of two tracers is fully implemented in our 
cornplex covariance matrix as off-diagonal terms. The result 
in iQ] would correspond to TZ/a ~ 1.8 at z < 1. Finally, 
any degeneracy with cosmological parameters in measuring 
the relativistic effects is largely eliminated, since the underly- 
ing matter distribution is cancelled iH. 

Primordial non-Gaussianity.— We extend our formalism 
to the primordial non-Gaussian signature in galaxy bias lfl9ll 
and introduce additional parameter /nl- Here we only con- 
sider the simplest local form of primordial non-Gaussianity 
to demonstrate how it can be implemented in the general rel- 
ativistic description, and ignore scale-independent and scale- 
dependent corrections (see, e.g., ll20ll ). 
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no. 2: Predicted constraints on the primordial non-Gaussianity /nl 
from galaxy power spectrum measurements. To facilitate the com- 
parison, the constraints on /nl are obtained by using the same sur- 
vey specifications as in 1/ ~ 50 (/i^^Gpc)'^ centered at jz = 1. 
Various curves show a /j^l with different priors on e and p. 

The primordial non-Gaussian signature in galaxy bias can 
be readily implemented in our full general relativistic descrip- 
tion with the Gaussian bias factor in Eq. (|2|i replaced by 

6^fe + 3/NL(6-l)^.^^g^, (6) 

or equivalently by adding additional term 3/nl(^ — 
l)5c^m{z) /T^{k, z) to the coefficient V in Eq. (|3]l, where 
T^f, is the transfer function of the curvature perturbation (see 

also Hill [H). 

In obtaining the constraint (Jf^^ on primordial non- 
Gaussianity, we set ctj. = c-p = 1 and marginalize over e 
and p. With the current uncertainties in e and p, the constraint 
CT/jjL (solid in Fig.|2]i is nearly identical to the unmarginalized 
constraint. The dashed curve shows that even with no pri- 
ors on e and p, <t/„l is not inflated except in the regime with 
CT/jjL ^ 2, because TZ and V are affected simultaneously by e 
and p but only V by /nl- Furthermore, the unique dependence 
of /nl on — 1 and Tip in Eq. ^ provides the multi-tracer 
method with more leverage to separate it from the general rel- 
ativistic effect. 

Finally, we allow e and p to vary as a function of mass 
with two logarithmic slope parameters and Up'. e = 



eol + ae ln(M/Mo), p = pol + ap ln(M/Mo) with cq = 3, 
Po = 0.4, ae = ap = 0, and Mq = IO^^H-'^Mq. The ef- 
fects of ae and ap (dotted in Fig. |2]i are sufficiently different 
from that of /nl, and (T/„l asymptotically reaches the floor 
set by the uncertainties in e and p. This demonstrates that 
the general relativistic effects in the galaxy power spectrum 
are not degenerate with the primordial non-Gaussian signa- 
ture. However, if /nl were to be constrained below unity, 
similar precision needs to be achieved in predicting TZ and V. 
This requirement can be relaxed by using independent mass 
estimates, as p drops out in Eq. ([3]l, further separating e from 
/nl- 

Discussion. — We have demonstrated that in a typical galaxy 
survey at low redshift general relativistic effects in the galaxy 
power spectrum can be measured at very high significance by 
using the multi-tracer method, providing a unique opportu- 
nity to test general relativity on horizon scales. We also show 
how the primordial non-Gaussian effect in galaxy bias can be 
implemented in the full general relativistic description, and 
quantitatively proved that the ability to detect primordial non- 
Gaussianity is little compromised by the presence of general 
relativistic effects. 

Proper application of our method to observations will re- 
quire a few key investigations beyond our simple treatment 
considered in this Letter. First, while we treated multiple 
galaxy samples as halos in multiple mass bins, it is shown 
11611 that a large scatter in the mass-observable relation, i.e., 
cinA/ = 0.5, would degrade the signal-to-noise ratio by less 
than a factor of two. Second, our prediction is based on 
the halo model description of the shot noise matrix, which is 
tested only for halos at M > 10^'^ Mq, and our prediction 
at M < 10^^/i~^A/o is an extrapolation. 
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